Negative radiation pressure exerted on kinks 
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The interaction of a kink and a monochromatic plane wave in one dimensional scalar field theories is studied. 
It is shown that in a large class of models the radiation pressure exerted on the kink is negative, i.e. the kink 
is pulled towards the source of the radiation. This effect has been observed by numerical simulations in the 
4> 4 model, and it is explained by a perturbative calculation assuming that the amplitude of the incoming wave 
is small. Quite importantly the effect is shown to be robust against small perturbations of the 4> 4 model. In 
the sine-Gordon (sG) model the time averaged radiation pressure acting on the kink turns out to be zero. The 
results of the perturbative computations in the sG model are shown to be in full agreement with an analytical 
solution corresponding to the superposition of a sG kink with a cnoidal wave. It is also demonstrated that the 
acceleration of the kink satisfies Newton's law. 

PACS numbers: ll.10.Lm, 11.27.+d 

I. INTRODUCTION 

It is by now universally accepted that spatially localized solutions of nonlinear equations (solitary waves, particle-like objects) 
are of great importance in several areas of physics, we refer to some of the recent monographs [0], |0], |0]. The study of localized 
solutions in various systems in a single spatial dimension (i.e. in 1 + 1 dimensional space-time) has proved to be quite fruitful in 
hydrodynamics, condensed matter physics and also for particle physics as a theoretical laboratory. A special class of non-linear 
equations admits "genuine" soliton solutions that retain their shape even after interactions, this feature distinguishes them from 
more generic particle-like objects. One can mention the sine-Gordon (sG) equation as a prototype example admitting "genuine" 
soliton solutions. It is of considerable interest both from a theoretical point of view and for physical applications to study the 
interactions of these spatially localized objects (not only of "genuine" solitons). There is a well developed framework for soliton 
perturbations in integrable and in near integrable systems 0], 10], where many powerful methods of integrable systems are 
available. In generic models, however, one has to resort to various perturbation techniques and in general also to numerical 
simulations. The motion of a localized object when subjected to a force has been shown to be governed by Newton's law for 
non-relativistic velocities |4]. It has turned out that in some cases, such as the sG model, deformation effects are also important 
and can lead to deviations from Newton's law 10], |0]. 

In this paper we study the motion of particle-like solutions in 1 + 1 dimensional scalar field theories (these objects are 
commonly referred to as "kinks") under the influence of an incident wave. In one of the most studied scalar models, in the "0 4 " 
theory, we have found that the kink starts to accelerate in the direction of the incoming wave. This effect has been first observed 
in Ref. in the present work we study it in much more detail both numerically and analytically. The unusual behaviour of the 
kink can be interpreted as being caused by a "negative radiation pressure" exerted on it. 

We have computed the force exerted on the kink by the radiation in a generic model in perturbation theory, assuming that the 
amplitude, A, of the incoming radiation is small. The leading order force exerted on a kink by a wave coming from the right 
(i.e. from x = oo) is then given as F = —A 2 q 2 |R| 2 , where q is the wave number of the incident wave, and R is the reflection 
coefficient. This means that the kink is pushed back by the radiation as expected. Now in a class of theories containing among 
others the <f) A and the sG models, the leading order reflection coefficient is zero, R = 0. In such cases the force is determined 
by higher order terms, in fact in the next order the force is F ~ 0(A 4 ). The basic physical effect responsible for the negative 
radiation pressure can then be understood as follows. For small enough amplitudes of the incoming radiation, the kink of the 
(f> 4 theory is transparent to the waves to first order in A. This is due to the reflectionless nature of the effective potential in the 
kink background. Because of the nonlinearities, during the interaction higher frequency waves are also generated. Some of the 
energy of the incoming wave is transferred to higher frequency (mostly double frequency) waves. These higher frequency waves 
carry more momentum than the incoming one. In the case of the </> 4 -type models the amplitude of the transmitted waves with 
double frequency is larger than that of the reflected ones. This way a momentum surplus is created behind the kink, thus pushing 
it forward. For large enough amplitudes the radiation pressure becomes positive. 

Comparing the perturbative result with numerical simulations, reasonably good agreement has been found for both the non- 
relativistic acceleration of the kink, and for the force acting on it. Also for small enough amplitudes of the incident wave the 
kink starts to accelerate according to Newton's law. 

There is no such negative radiation pressure exerted on the kink in the sG model, although the sG kink is also transparent 
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to the waves. As it turns out the sG kink is transparent to all orders in A (this special feature of the sG model is due to its 
integrability). In fact there is an analytical solution found in Ref. II 10TI corresponding to the nonlinear superposition of a kink and 
a traveling wave, which corresponds precisely to the problem we study. We have compared our perturbative results for the sG 
kink with this analytical solution and the perfect agreement found serves as a good test of our approach. 

In our context the negative radiation pressure appears because the kink is transparent to small amplitude incoming waves to 
first order in perturbation theory. Since this happens only in rather special cases, the structural stability of the effect should be 
addressed. We have demonstrated the robustness of the negative radiation pressure on the example of a generic perturbation of 
the </> 4 model. For an arbitrarily small perturbation of the cf) A theory, to first order in the amplitude of the incoming wave the 
radiation pressure exerted on the kink becomes positive, just as expected in a generic theory. The higher order contributions 
tend to compensate this however, and we have found that even for non-infinitesimal perturbations of the <fi 4 model there exists 
a critical amplitude of the incoming wave, above which the radiation pressure becomes negative again. This robustness of the 
effect makes it worthwhile for further studies. 

The organization of the paper is as follows. In Section II we introduce the models and present the first order calculation of the 
radiation pressure exerted by an incoming wave on the kink. In Section III we discuss the higher order perturbative calculations 
of the force and of the kink's acceleration. In Section IV these results are applied to the 4 and to the sG model, where an 
analytic formula for the force exerted on the kink is derived, and we also elucidate the physical reasons for the negative radiation 
pressure. In Section V the basic setup for the numerical simulations used is given and the analytical results are compared to 
those obtained by the numerical simulations. In Section VI the structural stability of the negative radiation pressure with respect 
to perturbations of the Lagrangian is studied. 

Most of the computational details are relegated to three Appendices. In Appendix A some details of the higher order pertur- 
bative computations are given, In Appendix B we perform a suitable expansion of the analytical solution of the sG kink on a 
cnoidal wave pertinent to our problem. 



II. FORMULATION OF THE PROBLEM 



A. The models considered 

We consider the following class of scalar theories in 1+1 space-time dimensions specified by the Lagrangian 
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c = ^d^d^~u{^), CD 

where the self-interaction potential, U (4>), is assumed to possess at least two degenerate minima (vacua), denoted by U (4> vac )). 
The equation of motion obtained from the Lagrangian ([TJ can be written as 

0-ct>" + U'(4>) = O, (2) 

The class of theories (Q]i admits static, finite energy solutions of the field equations (|2j. These solutions interpolate between 
two vacua of U((f>) and are commonly referred to as a kink or antikink. A kink (or antikink), (j) s (x), is given explicitly by the 
following formula 



x — xq = ± / — . . (3) 

J^xo) x/2C/(</)) 

In this paper we shall mostly work with the "0 4 " where U(<f>) = (<fr 2 — l) 2 /2, resp. with the "sine-Gordon" (sG) where 
U((f>) = 1 — cos0 models. Note that in the case of the </> 4 model, <fr vac = </>(±oo) = ±1, and 4> vac = </>(±oo) = 27r, for the 
sG equation. 

The corresponding kink solutions are well known: 

4> s = tanh x , in </> 4 , resp. 4> s = 4 arctan exp(— x) , in the sG model . (4) 

These kinks are well (exponentially) localized, and their position is conveniently defined by the zero of <fi s (x), which also 
corresponds to the maximum of their energy density. 



1 our conventions are: g M „ = diag(l, —1), dof = dtf = /, dif = 8 x f = /' 
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The energy-momentum tensor of the scalar field theory (Q~|i is given as: 

T^v = d^d v (f) - g^C , (5) 
whose components are spelled out explicitly for later convenience: 

Too = £ = \tf + \^ + U(<f>), (6) 
Toi = T w = -P = <l>'4>, (7) 
Tn = l^ + l^-U^), (8) 

where £ and V are the energy and momentum densities, respectively. The energy-momentum conservation laws d^T^ = can 
be written as: 



d t £ = d x , (9) 

d t V = - X -d x (ft + 0' 2 - 2U {</>)) . (10) 

B. Interaction of the kink with radiation 

The physical problem we wish to study is the interaction of a kink with an incoming (scalar) radiation. We assume that at least 
in the case when the radiation can be considered as a small perturbation, it is a reasonable approximation to treat the kink as a 
particle accelerating under the force coming from the radiation pressure exerted by the radiation. More precisely we consider 
the problem that a monochromatic wave, x), coming from the right (from x = 00) is incident on an initially static kink at 
x = 0. The incident wave x) itself reduces asymptotically to a plane wave, i.e. 

as — > 00) -> ARe{e* iuJt+qx) } , (11) 

where A is the asymptotic amplitude. The dispersion relation between the frequency, w, and the wave number, q, is easily read 
off from Eq. 

~u; 2 + q 2 + U"(^ vac ) = 0. (12) 

In order to allow for a perturbative solution of the time dependent problem, we shall assume that the amplitude of the incoming 
wave, A, is sufficiently small and expand the solution of the nonlinear wave equation (O in power series of A: 

= S + £ = tj> 3 + A^ + A 2 ^ + ■■■ , (13) 
To first order in A, the "radiation", i^ 1 ) satisfies a linear wave equation in the background of the kink: 

^)-^)" + u"(M^m (1) = l {1) + U {1) = o, (i4) 

where U"((f> s (x)) corresponds to the potential of the Schrodinger-type operator L = —d 2 /dx 2 + U" (<p s (x)). This potential can 
be written explicitly for the "(j) 4 " and for the sG models as 
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U"((j)f) = 4 — in the 4 resp. U"((j) s s G ) = 1 — in the sG model. (15) 

cosh x cosh x 



Introducing 



= \ {e wt v q (x) + e-*-V«(aj)) ■= Re^Me^} , (16) 



equation ( fT4l can be separated, where the r] q are eigenfunctions of the Schrodinger operator L, 

dx 2 



Ll h ■= [ -4-o + U"{4> s {x)) ) Vq (x) = u;%(x) . 
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As it is well known from elementary quantum mechanics, for potentials tending to zero for |ar| — » oo, the asymptotic forms of 
the scattering eigenfunctions are: 

^(a-^+oo) = < • n. iqx . (18) 

77,(2:^-00) = Te^, (19) 

where R and T are the reflection and transition coefficients. The first order solution, corresponds to the incoming radiation 
field, which reduces asymptotically to a monochromatic plane wave coming from the right to the kink. Let us remind the reader 
at this point, that the kink in any translationally invariant theory possess a discrete eigenfunction with eigenvalue uj = (a zero 
mode), which is called the translational mode. In some cases the Schrodinger operator L also possess other discrete eigenstates 
for lu > 0. There is one such a discrete (or internal) mode in the "(/> 4 " model for ujd = v3, and there is none in sG model. 

From the energy conservation law, Eq. (O it is easily seen that to first order the change of the total energy of the kink+radiation 
system in a box of size 2L, is 

L 

d t E = J dxdtS = A 2 e (1) 'e (1) \- L • (20) 

-L 

Averaging in time over a period, T = 2tt/uj, one finds 

(d t E\ L _ L ) T = A 2 qu J (l-\T\ 2 ~\R,\ 2 )/2, 

where (F)t denotes the average of the quantity F in time (over a period T). 

Assuming that the size of the box is sufficiently large and that the radiation of the kink itself can be neglected, the energy 
contained in the box is conserved, and then it follows that 

|R| 2 + |T| 2 = 1. (21) 

The rate of change of the total momentum in the box [-L, L] can be identified with the total force exerted on the the system of 
kink+radiation in its inside. To linear order in perturbation theory from Eq. (TlOb this force is found to be 

d t P\ L _ L ) =F^ = lA 2 q 2 (-l-\R\ 2 + \T\ 2 ) =-A 2 q 2 \R\ 2 . (22) 



- l /t ' 2" 

We now show that for not too long times the kink obeys Newton's law under the action of the force given by d22l . In 
order to define the acceleration of the kink we make the (usual) assumption that for small enough velocities one can neglect 
deformation and other effects and approximate the moving kink by cj>(x, t) = 4> s (x — X(t)). This approximation corresponds 
to introducing simply a collective coordinate for the position of the kink. Therefore for small enough displacements, <f>{x, t) m 
4> s (x) — X(t)(j>' s (x). In fact r)t ■= 4>' s (x) is nothing but the translational zero mode of the kink, which is orthogonal to all other 
linearized (internal and radiation) modes of the operator L. The leading order displacement of the kink is then easily obtained: 

X(t) = -j^, (23) 

where (f\g) denotes the natural Hilbert space scalar product. The acceleration, a to leading order can then be calculated as 
follows 

o(w , = jt (t) = Jnm = JvS^t&i) (24) 

where ^ n > is the lowest order approximation i.e. the smallest n for which the acceleration is non zero. To compute higher order 
corrections to the acceleration is nontrivial, since in our approximation the time dependence of the kink has been encoded in 
the single collective coordinate, X(t), whereas one has to take into account distortion, radiation, etc effects. In the generic case 
discussed above n = 2. Therefore to lowest order in PT the force ( f22l is quadratic in the amplitude, A, and clearly according 
to our definition (124-b the acceleration of the kink is also 0(A 2 ). Then the second order perturbative solution, ^ 2 >, is needed. A 
not too difficult computation (see Eqs. (IA7K (IA1 lb in Appendix A) yields: 

m s a^ = -A 2 ( Vt \^ 2) ) = ^(tf'"(&)»W7-,|Tfc) = I*!' ■ (25) 
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where the relation (i]t\i]t) = ms has been used. The method to project unto the translational mode to compute the kink's 
acceleration was also used in Ref. [&] where the dynamics of a kink in the <fi 4 model with a perturbed potential was studied. 
One sees that to leading order in PT the acceleration of the kink, a*- 2 -*, is indeed given by Newton's law, = m s a^ 2 \ In 
particular one can see that the kink is pushed back under the action of the force coming from the radiation pressure as expected. 
Therefore it is consistent to identify the time averaged momentum flow in the box Eq. d22b . with the total force acting on the kink 
to leading order. When the radiation field can be treated as a small perturbation the effect of the momentum flow on the radiation 
field itself can be neglected as a first approximation. One would then expect that for small enough amplitudes of the incident 
radiation, it is a reasonably good approximation to the solution of equation (0, that an initially static kink starts to accelerate as 
a non-relativistic particle of mass m s . It is natural to expect that other effects, such as the radiation by the kink, its distortion, 
etc. show up only in higher orders. We remind the reader that identifying the force exerted by the radiation on the kink with 
the total momentum flow in the box is valid only to leading order (and after time averaging), and also for relatively short time 
intervals. The force ( f22b is quadratic both in the amplitude of the incoming wave and in the reflection coefficient, in complete 
analogy to the well known radiation pressure in classical electrodynamics. 

It is quite illuminating to compare the prediction for the acceleration of the kink (f2]) by solving numerically the nonlinear wave 
equation (O for a few common one dimensional field theories. On FigureQ]we have depicted the trajectories of the (topological) 
zeros of various kink solutions under the influence of an incoming radiation from +00 in three models: in the sine-Gordon, cf> 4 
and in the <ft 8 model where U = \(<fi 2 — l) 4 . It is somewhat surprising that the kink has been pushed by the radiation pressure 
only in the very last example, in agreement with Eq. d22l i . As it can be seen on Fig. Q]the time average of the acceleration of the 
kink in the sG model is zero, the kink is steadily oscillating around its initial position. Most remarkably the kink in the <f> 4 model 
accelerates towards the source of radiation, and it is this interesting effect that we interpret as negative radiation pressure. Taking 
into account other collective coordinates such as the shape mode, would not substantially influence our main results, therefore 
we have chosen to ignore them. Note that the acceleration of the kink in the <fi 8 model is noticeably larger as compared to the <fi 4 
one, for the same amplitude of the incoming wave (A = 0.14). In fact, while according to Eq. (|22| | the acceleration of the kink 
is quadratic in the amplitude A, this is only true for the </> 8 model. The acceleration of the <\> 4 kink turns out to be proportional to 
A , Motivated by the unexpected results shown on Figure[T]we shall investigate the interaction of a kink with radiation in more 
detail in order to give an explanation of the negative radiation pressure. It is immediately clear that the derivation of the force 
obtained in Eq. (f22b is only valid if its leading contribution comes from the linear approximation. In the case when |R| <C 1 the 
validity of this assumption is questionable since then the higher order terms may contribute in an important way, therefore the 
first order result in Eq. d22l i is not necessarily correct. As a matter of fact it is rather well known that in both the <fi 4 and in the 
sG models that the potentials ( TT~5T > are reflectionless, i.e. R = 0, for all frequencies. Therefore in such models the dynamics of 
the kinks is determined by higher order terms. 



III. COMPUTATION OF THE FORCE ON THE KINK UP TO 0(A 4 ) 



In this Section we shall outline the computation of the force acting on the kink as well as its acceleration to higher orders of 
perturbation theory. As we shall show below the next nontrivial correction to the force turns out to be of 0(A 4 ). Fortunately 
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one does not have to go up to 4th order in perturbation theory. 

By expanding the equation of motion of the kink (0 in power series of the amplitude of the incoming wave, A, corresponding 
to ( fT3l ), the nth order solution of the equation of motion is determined by the inhomogeneous linear equations 



00 



00 



/ 



00 



(26) 



where the source terms /W can be calculated from the lower order terms in the perturbation series, see AppendixlAlfor details. 
To define the solution of the perturbative equations d26b uniquely, we impose in each order that there are no incoming waves 
from the left (x < 0), as boundary conditions. These correspond to the physical problem of an incoming wave from the r.h.s. of 
the kink. We remark that to avoid resonances in higher orders, the frequency, u), has also to be expanded as 



(27) 



The total force acting on the system "kink + radiation" inside the segment [-L, L], can be calculated to higher orders in A 
from the energy and momentum conservation laws similarly to the leading order computation. We compute the rate of change 
of the total momentum inside the box [-L, L] to obtain the total force acting on the system, which after averaging in time can 
be consistently identified to leading order with the force acting on the kink just as in the calculation in the previous section. 

Using <f)(x, t) = (j) s + £ in Eq. dTob and integrating over the interval [-L, L], one finds 



d t P 



1 



(28) 



where we have omitted terms which are exponentially small for L \. Choosing L to be sufficiently large, one can use the 
asymptotic form of the higher order solutions to find the rate of change of the averaged momentum inside the segment [—L,L] 
just as for the lowest order case. Assuming that dtP = to 0(A 2 ) in Eq. d28l l (reflectionless case) one finds that all terms of 
order 0(A 3 ) also drop out, i.e. there is no momentum flow to this order into the segment (for more details see Appendix A). 
Therefore the first non-trivial contribution to Eq. (l28T l is of order 0(A 4 ), implying that one would also need £' 3 ', i.e. one should 
compute up to 3rd order in perturbation theory. This complicates considerably the problem, even if only the asymptotic forms 
of the are needed. For some computational details of the higher order calculations we refer to Appendix A. Remarkably 
though one can in fact eliminate the contributions coming from the 3rd order terms from the momentum balance in the segment 
[-L, L] by exploiting the law of energy conservation. As it will be shown below all the information needed to calculate the force 
acting on the kink is actually encoded in the asymptotic form of the time-dependent part of £^ 2 ) . 

Denoting by (ffl the m-th coefficient in the Fourier expansion of £ in the 0(A n ) order, the 0(A 2 ) order solutions can be 
written as: 



(2) 



J2lLOt 



t(2) 
S+2 



t(2) 
SO 



-2iujt 



t(2) 
S-2 



(29) 



(2)* 



We also note that for reflectionless potentials £ +2 
term, £( 2 ), has the form 



(2) 

£_2- The asymptotic form of the time-dependent part of the 2nd order 



4+2 ( x ->■ ±°° 



U'"(ct> vac ) 



rf+ q + a 2 2,±k(q)VTk , 



where 



24E/"(</w) '+« 

i r°° 

a 2 2A<l) = -^ J_ dx> Vk if q U"'(<Ps) 



(30) 



(31) 



k = \J AbS 1 — U"(4>vac) is the wave number corresponding to 2u>, and W = ilkV-k ~ VkV-k = —2ik is the Wronskian. In a 
way a 2 2,k(q) encodes the reflection and transition coefficients due to the non-linear effects. The details of the computation of 
Eq. d28| ) up to 0(A 4 ) can be found in Appendix A, leading to the result 



( dt P) T = F< 4 > = -A 4 [2k 2 (\a 



22, + fc I 



*22, 



s |) -2<7 2 Re(a3i,_ 9 ) 



(32) 



where the pertinent contribution from the 3rd order terms is encoded in a single coefficient 031,-5, defined in Eq. dA25l ). In Eq. 
(T32t we have identified the time average of the overall momentum, (dtP)r, flowing into the segment [-L, L], with the force, F, 
exerted by the incoming radiation on the kink. 

In the case of reflectionless potentials the first non-vanishing contribution to the acceleration comes from the 0(A A ) terms. 
Considering for simplicity such theories where the kink is spatially antisymmetric (this includes both "0 4 " and the sG models) 
a straightforward computation yields: 



(33) 
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i.e. it is sufficient to compute the second and third order solutions (see Appendix [A] for more details). By a direct computation 
we have checked that with the definition of the force in Eq. (l32l acting on the kink within our perturbative framework, Newton's 
law 

F (4 » = m s a^ , (34) 

still holds, at least up to the fourth order 0(A 4 ). In our view this result lends strong support to identify the time average of the 
momentum flow in the segment [-L, L] with the force acting on the kink to leading order in PT. 

We give here a derivation of Newton's law, which also indicates the limits of its validity. Using the following simple identity 

= dt{Hf - &') = -d t (V + &') , (35) 
by integrating over a segment [— L, L] one easily obtains the relation 

L L 



J dx'irit =-Pl- J dxdt{iC) , Pl = P 

-L -L 



-L 



(36) 



where Pl denotes the momentum change inside the segment [— L,L]. We shall approximate the lhs of Eq. d36t simply by 
(tlVt) : = —m s a since the difference between them is exponentially small in L. After averaging in time we obtain 

L 

Fl =m s a- J dx (($(&')) T , (37) 

-L 

where F L denotes the total force acting on the box. Now in perturbation theory the solution £ can be decomposed as 

£ = C P (t, x) + £ (x) - ^at 2 £ t {x, t) (38) 

where £ p (t, x) is periodic in time and it is at least of order O(A), £t corresponds to the "accelerating part" (with initially constant 
acceleration), and £o(x) is the time independent piece. This holds to order 0(A 2 ) when ^ 0, resp. to order 0(A 4 ) when 
o« = 0: 

£ P = A^ + A 2 e (2) + 0(A 3 ) £ (x)=A 2 tf\x) + O(A 4 ) 6 - A n Vt (x) + 0(A n + 2 ) , (39) 

where n = 2 if ^ 0, and n = 4 if a*- 2 ' = 0. The leading order correction is oc aA n+1 , which depends on time averaging 
and indicates the limits of the validity of our simple minded approach. 

Let us note here that the mass of the kink gets renormalized due to its interaction with the radiation field. A standard calculation 
for the first correction to the kink mass is given as m* = m s + A 2 5m^ ■ • • , i.e. to lowest nontrivial order 5m is proportional 
to A 2 , therefore to 0(A i ) it does not show up in Eq. fl34l i. Nevertheless the numerical simulations (see SectionO indicate that 
the kink has an effective mass, whose value is quite close to the renormalized mass. 

The computation of the rate of change of the energy inside the segment is completely analogous to the previous momentum 
balance calculation, we find: 

(d t E) = -A 4 [Auk (|a 2 2:fc | + |a 2 2 ,_ fc |) + 2c^Re(a 3 i,- 9 )] ■ (40) 

Assuming that after averaging, at least for some initial time the kink can be considered as a rigidly accelerating particle, if it was 
initially at rest, i.e. v(t = 0) = then obviously 

(d t E}\ t = = mvv\ t = = . (41) 

This equation together with eq. d40b can be now used to eliminate the coefficient 0:31,-5 from Eq. d32l , and then one obtains a 
remarkably simple formula determining the force acting on the kink: 

F<$ = 2A 4 k [(k - 2q)\a 2 22 ^ k \ -(k + 2q)\ a 2 ^ k \} . (42) 

By a direct computation of the energy we have verified that (dtE)\ t= o is indeed zero up to 0(A 6 ), which shows the validity of 
Eq. (SB. 

To conclude this Section, we have calculated the force exerted by an incoming wave on the kink to the first non-trivial order in 
perturbation theory in the class of models where the linearization around the kink yields a reflectionless potential. It is important 
to emphasize that for the class of models where the effective potential is reflectionless, the force exerted on the kink turns out 
to be proportional to F ~ 0{A A ). This is to be contrasted to more generic models where the effective potential is reflective, in 
which caseF~ C(^l 2 ). 
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IV. NEGATIVE RADIATION PRESSURE IN THE <j> A MODEL 



In this Section we apply the previously obtained general results to compute explicitly the force exerted on the kink by an 
incoming wave in the <f> 4 and in the SG models. 

In the cj) 4 model the full nonlinear equation for the "radiation" £ is: 



£ + U + 60 s <e 2 + 2£ 3 = . 
The first order solution £W in Eq. ( TToT l can be explicitly given both in the <fi 4 and in the sG models: 



and 



3 tanh x — 1 — q 2 — 3iq tanh x 

V(g 2 + l)(<z 2 + 4) 

sG *9 — tanhx 



e %qx , where q 2 + 4 = lu 2 , 



e tqx , where q + 1 = lj- . 



(43) 



(44) 



(45) 



The 2nd order solution (0(A 2 )) given by Eq. d29b contains the zero frequency term, Q , which in general depends both on t 
and x. In the present case it is consistent to assume that £q 2 ' is time-independent (see in Appendix lAl for a proof), moreover 
it can be written explicitly. The 2nd order "transition and reflection" coefficients, a22,k(q) in Eq. ( f3TT > can also be calculated 
analytically, the result is: 



3 q 2 +4 jq 2 + 4 

0122 Ml) = "o 71 " 



1 



2 q 2 + l\ k 2 + l ksinh (i3±]l w ) ' 



where q 2 



(46) 



4, k 2 = 4(uj 2 — 1). Using Eq. d46b it is now easy to find the averaged force d42b exerted on the <p 4 kink by an 



incident wave of frequency to, we obtain 



^(4) 



9ir 2 A 4 (j 6 



fc(4w 2 -3)(w 2 -3) 5 



U>- 



<jj+ 



sinh 7tcj_ sinh 7rwj 



with 



:= V^ 2 



1 ± vV 



(47) 



Introducing = /A 4 , 
2while/( 4 )(w; 



the behaviour of the function f^ 4 ' [to) for uj — > 2 (i.e. small values of q) is given as /^ 4 -* (u> — * 2) f=a 
^3/4 for large values of w. Quite interestingly the force acting on the kink is positive, therefore 
This is the effect we refer to as negative radiation pressure. The origin of the 



0.3749V^T 

it accelerates towards the source of radiation, 
negative radiation pressure can be understood by noticing that for all frequencies, u>, the amplitude of the (non-linearly) reflected 
wave |a22,+fc| is smaller than the amplitude of the transient wave |a22,-fc|. i- e - a 22 k ^ a 22 -k- ^ n ^ rst order perturbation 
theory, such an effect would not be possible because of the identity ( T2"TT i expressing energy conservation at the linear level. 

The surprising effect of negative radiation pressure on the kink exists only because of the presence of nonlinearities. In the 
linear approximation the kink is transparent to the incident wave, therefore it does not accelerate. Because of the nonlinear terms, 
part of the energy of the incoming wave is transformed into a wave whose frequency is twice of the original one. This double 
frequency wave has a larger ratio of momentum to energy density than the incident wave with smaller frequency, hence it carries 
more momentum than the originally incident one. This way a surplus of momentum is created behind the kink, which pushes 
it towards the direction of the incoming wave. The above is of course only an intuitive explanation of the negative radiation 
pressure on the kink. The effect of negative radiation pressure has been clearly observed in our numerical simulations of the <fi 4 
model (see the Section "Numerical simulation"). 

Next we compute the acceleration of the <f> 4 kink using Eq. (|33T l derived in the previous Section. We need to compute in 
fact the projection of the second and third order solutions on the translational mode, ij t . We have computed them numerically 
by two different methods, using the integral representation based on the explicitly known Green's function and also by direct 
numerical integration of the corresponding equations d26l i. We show separately the three projections in ( l33l on Figure[2](divided 
by (j] t \r]t) = 4/3). It is worthwhile to point out that all the three projections are positive. 

Let us now turn to the sine-Gordon kink, and evaluate the force of order 0(A 4 ). In the case of the sine-Gordon model our 
numerical results did not show any net radiation pressure, the sG kink was oscillating around its initial position (c.f. Fi^JJ). 
Interestingly in the sG model the second order waveform can be calculated in closed form: 
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FIG. 2: Three projections obtained from the second and third order (needed for!33b. 



from which one can immediately see that the coefficients 022, ±k determined by the asymptotic behaviour of £( 2 ) in Eq. (f30b are 
zero. This implies that the radiation exerts no force at all on sG kink at least up to this order. Clearly this interesting fact should 
be related to the special feature of the sG model, namely its integrability. As a matter of fact there is a remarkable analytic 
solution corresponding to the nonlinear superposition of a kink with a travelling (cnoidal) wave in the sG model obtained by 
H. J. Shin ITotl using the Darboux transformation method. In Appendix Icl we give a short review of Shin's solution which is 
somewhat complicated, and we demonstrate that one can expand it in a parameter which can be identified with the asymptotic 
amplitude of the incoming wave. This way we have verified that both the first and the second order solutions obtained by our 
perturbative calculations agree perfectly with the Taylor expansion of the analytical solution. This comparison has also served 
as a test of the validity of our perturbative method. 



V. NUMERICAL SIMULATION 



In the present Section we outline the numerical method used to solve the nonlinear wave equation Eq. (O describing the 
interaction of a kink with an incident wave and present the results in form of figures and tables. 

We have discretized equation in the spatial variable, x as cj>(nh,t) := (f> n {t). The second derivative was approximated 
using the following five point scheme: 



12h 2 



( 



30^„ + 160. 



'ra+l 



0„ +2 ) + 0(h 4 ) . 



This way Eq. (0 reduced to a system of ordinary differential equations: 

4> n = dk - u'{<j> n ) . 



(49) 



(50) 



We have simply put this coupled infinite system into a finite box of size 2L, which was then solved using a standard fourth order 
Runge-Kutta method. 

Our initial conditions have been chosen to correspond to a kink together with a first order travelling wave: 



(j)(x,t = 0) = <j) s (x) + -Arjg (x) +c.c. , 

4>(x,t = 0) = -iuiAr]q(x) + c.c. , 

and we have fixed the boundary values of <p(x, t) at x = ±L as 

<f>(x = ±L,t) = ±1 . 



(51) 
(52) 

(53) 



The evolution time of the system was restricted to be smaller than L, to avoid the unphysical influence of the reflected waves from 
the artificial boundaries at x = ±L on the kink's motion. The position of a static kink can be quite unambiguously identified 
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FIG. 3: The position of the zero of <j> as a function of time for A — 0.12, to = 3.0. 



Amplitude A 


Fitted acceleration 


Theoretical value 


0.10 


0.0000482 


0.00004354 


0.12 


0.0000977 


0.00009029 


0.16 


0.0003188 


0.00028537 


0.18 


0.0005296 


0.00045710 


0.20 


0.0008325 


0.00069670 


0.22 


0.0012541 


0.00102005 


0.24 


0.0017300 


0.00144469 


0.26 


0.0014806 


0.00198986 



TABLE I: Fitted and theoretical values of the acceleration for u = 3.0 

by the location of its zero which coincides with the maximum of its energy density. It is less clear how to define the position 
of an interacting kink. In our case one has to separate first the field of the kink from that of the radiation, which can already be 
problematic and the position of the kink is not very well defined, in general the maximum of the energy density and the zero of 
(j)(x, t) do not coincide. For small enough amplitudes the kink is only slightly perturbed and therefore its topological zero is still 
a rather satisfactory definition as the position of the kink and we have used this definition in our work. 

We have plotted the position of the zero of 4>(x, t) as a function of time for the frequency u> = 3.0 and for the amplitude of the 
wave A = 0.12 on Figure 3. On this figure one can clearly see that the trajectory of the zero of <fi(x, t) is quite close to a parabola, 
corresponding to the expected non-relativistic acceleration of the kink. A numerical fit confirms that a parabola of the form at 2 /2 
to the trajectory is good a approximation indeed. The fitted value of the acceleration in this case was a num = 9.72 • 10" 5 which 
is not that far from the result of our analytical calculations in Eq. d33l giving a t hoor = 9.02 • 10 -5 . Taking into account that in 
the analytical calculation only the leading terms have been used, this agreement appears to us satisfying. Next we have checked 
if the measured acceleration is indeed 0(A 4 ) as predicted by the leading order perturbative result ( [33) , On Figure|4]we have 
plotted the fitted acceleration for uj = 3.0 divided by A 4 for the amplitudes of incoming wave varying between 0.1 < A < 0.3. 
As one can see for 0.1 < A < 0.22 the curve is close to being fiat implying that the dominant term is indeed proportional to 
A 4 . This proportionality breaks down when the value of the amplitude increases to about A w 0.24. For A = 0.25 even the 
sign of the fitted acceleration changes. In Table I we have compared the numerically obtained values of the acceleration to the 
theoretical ones for a range of amplitudes. 



From Table I one can see that the agreement between the calculated and the fitted values of the acceleration is reasonably good 
up to values of A < 0.22. These results confirm that for amplitudes of the incoming wave in the range 0.1 < A < 0.22 the kink 
accelerates non-relativistically, and also that its acceleration scales as A 4 . Next we exhibit the numerically obtained acceleration 
in function of the frequency uj on Figure[5] together with the theoretical curve, and some results are given in Table II. 
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FIG. 4: Fitted acceleration divided by A 4 for lo = 3.0 as a function of A. 
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FIG. 5: Fitted acceleration divided by A 4 for ^4 = 0.16 (solid line) and the acceleration calculated analytically (dashed line). 



The first thing one might notice on Figure [5] is the presence of three resonance-like structures completely absent from the 
theoretical curve which is a monotonously increasing function of lo. The largest resonance is not very far from 2ujd ~ 3.46 
which indicates that it is likely to be related to the coupling between internal (or shape) mode of the kink and radiation. A 
plausible explanation of the important change in the acceleration at frequencies when the shape mode couples strongly to the 
incoming wave is the following. At such a "resonance" frequency the shape mode accumulates a substantial amount of energy, 
which is then radiated symmetrically in both directions. Far from the kink this produces the same effect as a reflected wave, thus, 
at such resonant frequencies the kink is not transparent. On Figure [6] the path of the kink is plotted for a value of the frequency 
near the resonance. As one can see the motion of the kink is somewhat irregular there. All in all the discrepancy between 
the results of the perturbative computations to leading nontrivial order and those of the numerical simulation does not exceed 
10% for a large range of the frequencies with the exception of three resonances. It seems to us that in view this agreement is 
satisfactory in view of the approximations used. 

As it can be seen from Table U the numerically found acceleration is systematically larger than the leading order theoretical 
one. Clearly higher order effects could play a role here, and the simplest one to be taken into account is the renormalization of 
the kink mass due to the radiation field. The lowest order 0(A 2 ) contribution to the mass, Sm^ is negative for all frequencies, 
i.e. the effective mass to* = to + A 2 8m^ < m which goes into the right direction. It is quite difficult to obtain a sufficiently 
precise numerical value for the effective mass, nevertheless we have obtained some indicative results. The numerically computed 
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Frequency u 


Fitted acceleration aA 4 


Theoretical value 


2.50 


0.4285 


0.339743 


2.70 


0.4472 


0.388489 


2.90 


0.4397 


0.421959 


3.10 


0.4865 


0.446130 


3.30 


0.3510 


0.464260 


3.50 


0.4761 


0.478267 


3.70 


0.3310 


0.489349 


3.90 


0.4322 


0.498290 


4.30 


0.4672 


0.505624 



TABLE II: Fitted acceleration divided by A A for A = 0.16 
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FIG. 6: The position of the zero of <f> as a function of time for u — 3.38 and A — 0.14. 

momentum balance (force) inside a box [-20,20] is presented on Figure [8] Dividing this force by the measured acceleration 
we obtain the effective mass, to*. The mass obtained in this case is m* um = 0.954to s . An analytical calculation yields 
— SmS 2 ' = 3A 2 (uj 2 — 2)/(uj 2 — 3), which gives m* heor = 0.962m s , so there is a reasonable degree of agreement between the 
two. In conclusion the numerical results show that our perturbative calculations are quite reliable for amplitudes A < 0.2 and 
for frequencies far from resonance points. Finally on Figure [7] the behaviour of the acceleration of the <p 4 kink for a range of 
amplitudes and frequencies of the incoming wave is depicted. 

VI. STABILITY OF THE EFFECT UNDER PERTURBATIONS 

In this Section we shall demonstrate that the effect of negative radiation pressure has a certain degree of robustness with 
respect to perturbations of the original </> 4 model. This fact makes the effect, which is in itself interesting, much more relevant 
for physical applications. At first sight it is not so obvious that this effect could survive a small perturbation of the model at all, 
since a generic perturbation, no matter how small it be, destroys the reflectionless nature of the potential in Eq. < fl~4T >. This way a 
first order perturbative contribution (in the amplitude A) is generated. Therefore the leading term in the expression for the force 
changes under the influence of a generic perturbation from being of order 0{A 4 ) in Eq. (l42l to 0(A 2 ) as in Eq. (l22t . We shall 
show that although a small, but generic perturbation of the <fi 4 model changes the leading term for the force to being of order 
0(A 2 ) indeed, for an important frequency range still the 0(A A ) will dominate if the amplitude A > A crlt . We shall consider a 
concrete example of perturbation which illustrates that the critical amplitude A cr it can turn out to be small and that the negative 
radiation pressure stays practically unaffected. 

To start with we shall consider a generic perturbation of the field equation Eq. (O of the form: 

4>-<P" + U' ((/))+ e5U'(<P) =0, (54) 
where e is a small parameter and 8U'(<f>) is the perturbation. We look for the solution of the perturbed equation d54T i again 
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FIG. 7: Acceleration aA 4 as a function of amplitude and frequency of incoming wave. 
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FIG. 8: Momentum balance inside the segment [-20,20] for u) = 3.0 and A = 0.12. The slope of the fitted straight line (as a sort of averaging) 
is —1.24 • 10 _J , corresponding to an effective mass m* = 0.954m. 



expanded in a power series in the amplitude of the incoming wave: 

( f> = Mx)+M (1) (x,t) + -- - , (55) 

where <j) s (x) denotes the static kink solution of the perturbed field equation d54b and ^'(x, t) satisfies the field equation lin- 
earized around <f> s (x), 

£(D _ £(D" + [/"(0 S )£ (1) + eSU"^)^ = , (56) 
Now we also expand the solution of Eq. ( f54b in e determining the perturbation of the original theory: 

Mx) = 4 0) + «f>P +■■■, = + < (U) + • • • , (57) 

where (j)^ (x) is the static kink of the unperturbed equation 10, and £( 10 ' (x, t) is a solution of the linearization of the unperturbed 
field equation around <j)f^ (x), i.e. £( 10 ) + L^ 10 ' = 0. The equations for the first order corrections in the perturbative parameter 
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e are: 



£(11) _ £(11)" + v n Uo)^ ^(11) + p,,, ^(0)^ (1)£(1O) + ^„ ^(0) j £(10) = _ 



(58) 
(59) 



Let us recall that £ (10) (x, i) - 
following equation: 



1 e llJJt r]q + c.c, and look for the solution as £( n ) = i e 



zoji^(l) 
2 - f/ 9 



c.c, where rjq (x) satisfies the 



L — a; ) „W + 



(60) 



Since ( f60b is an inhomogeneous equation of the form of ( IA1 3b its general solution can be obtained from Eq. (IA14b . In order to 



obtain the force we need the reflection coefficient, R, therefore it is sufficient to compute the asymptotic behaviour of r/q for 
large \x\, which can be found to given as: 



rj^ (x — > +oo) 

Tj^{x — > — OO) 



(0) 
V-q 

W 

(o) 

w 



dx 



dx 



SU" 
5U" ( 



U" 



u' 



gW (0) 
PR 'l-q 



WW 



(61) 
(62) 



The expressions and e/J^P are the first non-trivial corrections to the reflection (R) and transition (T) coefficients. Recall 
that in the reflectionless case, such as in the 4 model R = and |T| = 1. From Eq. doTt it immediately follows that in the 
perturbed field equation, d54l i the kink is not transparent anymore. The leading contribution to the force is then determined by 
the first order term linear in the amplitude A. As found in Eq. ( f22b the dominant part of the force acting on the kink, F^ 2 \ is 
proportional to the square of the reflection coefficient, i.e. 



F 



(2) 



A 2 /^ = -q 2 \R 2 \A 2 « -q 2 e 2 \(3 ( ^\ 2 A 2 



(63) 



The above equation holds if \ef3^ | , |e/3jj | <C 1, which is true for sufficiently small values of e and for a certain range of q. The 
first perturbative correction in e to the force of order 0(A 4 ) in the amplitude, F' 4 ' = A 4 /' 4 ) (c.f. Eq. (l42l ) will be at least of 
order 0(A 4 e 2 ) therefore this term can be neglected in the following. The leading contribution to the force acting on the kink 
due to a perturbation of a theory with reflectionless potential is given as 



F = A 2 (-q 2 e 2 \(3^\ 2 + A 2 f^). 



(64) 



Assuming that f( 4 > > (i.e. that the radiation pressure is negative in the model) it follows from Eq. d64l ) that the amplitude of 
the incoming wave must be larger than a critical value, A > vl cr it for a fixed value of e to ensure F > 0, i.e. that the effect of 
the negative radiation pressure be present. The value of the critical amplitude is determined by the condition F^ + F 1 - 4 ^ = 0, 
leading to 



Ac. 



(65) 



Clearly the result in (1651 1 is meaningful in our perturbative framework, only if j4 cr i t <C 1. 

We now apply the above general results to the 4 theory (where U'(<f>) = 2(f)((f> 2 — 1)). We have chosen the following 
perturbation for 5U'(4>): 



5U'&) = </# 2 



1) 



Recall that c/)^ (x) is nothing but the static kink in the unperturbed 4 theory, so 

^(x) = tanhx. 

The first order correction in e to the static kink is determined by Eq. ( 1581 . which can be analytically solved: 

1 tanhx 



6 cosh 2 x 



(66) 



(67) 



(68) 
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One can also evaluate the integral d6TT l analytically, and the first order result in O(e) for the reflection coefficient is given as: 

15sinn7rg 

Evaluating the critical amplitude in Eq. ( l65l l for some values of the perturbation parameter e, one finds for e = 0.1 and q = 1.2 

(corresponding to uj f=s 2.33) A cr ;t ~ 0.021 and A cr it « 0.041 for e = 0.2. For increasing values of the frequency to, /3j^' — ► 
exponentially fast therefore A cr a becomes very small. 

In order to confirm the existence of the critical amplitude and compare its magnitude with the one found in Eq. d65l ). we 
have performed some numerical simulations on the perturbed <fi 4 model (1661) . It is not an easy numerical task to measure the 
critical amplitude since the acceleration tends to be very small for A ~ A crit and therefore long-time evolution with large spatial 
resolution is needed. The value of the perturbation parameter, e cannot be chosen to be too small either since then A cr i t becomes 
so small that we cannot measure it. Also the initial conditions used in the simulations of the perturbed (j> 4 theory gave some 
non-zero contribution to the initial velocity of the kink. The measured values of the accelerations were as small as 10~ 9 , 10 -8 . 
Fore = 0.1 andw = 2.33 we have found that A crit = 0.0168 ±0.0001, and A CTit = 0.0292 ±0.0002 for e = 0.2. The measured 
values for A cr it do not agree very precisely with the prediction of Eq. d65l ), the discrepancies being of about 20% resp. 30% for 
e = 0.1 resp. for e = 0.2. Nevertheless we consider the numerically found values to be consistent with Eq. d65l ) for the following 
reasons. A cr - lt /e still varies considerably (by ~ 13%) for the two considered values of e, indicating that higher order corrections 
are still important here. For the smaller value of e the discrepancy between the measured value and the result of Eq. d65l l is also 
smaller. In any case we have been able to demonstrate the existence of a critical amplitude in the perturbed <fi 4 model above 
which the radiation pressure becomes negative. The measured value of A CT n is consistent with the theoretical estimate in Eq. 
(|6B1 . Moreover the negative radiation pressure persists for rather large values of the perturbation parameter, e. For example in 
our numerical simulations for A = 0.1 and ui = 3.0 the radiation pressure became positive only for either e > 1.3 or e < —1.2. 
We have shown that for a rather large range in the magnitude of a generic perturbation of the <f> 4 model, the phenomenon of 
negative radiation pressure persists. 

Finally comin g to the perturbation of sG model, we have shown that to order 0(A 4 ) the force is zero. In fact using the analytic 
solution of Shin iflQIl we conclude that the force is zero to all orders. Therefore our derivation for the critical amplitude above 
which the force acting on the kink is still negative does not apply to the particular case of the sG model. We have observed 
that the kink only oscillates around its initial position, and the average of its velocity is zero. This implies, that even a small 
perturbation of the sG model may change this qualitative behaviour and determine the motion of the kink. 



VII. CONCLUSIONS 



We have studied the interaction of a kink in 1 + 1 dimensional scalar models with an incoming wave in perturbation theory. 
We have shown that in a certain class of theories (such as the <\> 4 model), the kink is pulled towards the direction of the incident 
radiation, instead of being pushed back. This interesting phenomenon constitutes an interesting example of negative radiation 
pressure, which in this case is due to the nonlinearities and to higher order effects. Comparing the results of the perturbative 
calculations to numerical simulations in various field theoretical models (mostly in </> 4 and in sG) a rather good agreement 
was found for not too large values of the amplitude of the incoming wave (up to A m 0.2). We have also addressed the 
important problem of structural stability of the negative radiation pressure with respect to generic perturbations of the theory. 
In models where the reflection coefficient is small we have established the existence of a critical amplitude above which the 
kink experiences negative radiation pressure. This is closely related to the robustness of the effect which has been demonstrated 
on the example of a perturbation of the <f> 4 model. We have found in fact that even for large perturbation there exists a critical 
amplitude of the incoming wave above which the radiation pressure becomes negative. 

In the sG model the radiation pressure turns out to be zero. In this model there is an analytical solution corresponding to the 
superposition of a kink and an incoming (cnoidal) wave ifioll . and we could confirm the correctness of our perturbative results 
by comparing them to the expansion of the exact solution. 

We have also shown that under the action of the averaged force exerted by the radiation the kinks accelerate in all these models 
according to Newton's law. 

It is clearly an important open question if the effect of negative radiation pressure is also present in other, in particular higher 
dimensional theories. Our preliminary results suggest that this phenomenon is also present at least in two other models of 
quite some physical interest: in the 2 + 1 dimensional complex (f> 4 theory (Goldstone's model) and in the Abelian Higgs model 
admitting vortices. This suggests that this effect might not be so rare as one could have expected at first sight. 

The phenomenon of negative radiation pressure is relevant not only for the interaction of a single kink with radiation but also 
for system of many kinks. One would expect for example that a many-kink system where the radiation pressure is negative 
becomes unstable with respect to collapse. Since the interaction of well localized kinks is rather weak (e.g. for two <fi 4 kinks at 
a distance L the force between them is F ~ e~ aL ), the radiation pressure can play an important role in many kink systems. 

It would clearly be interesting to investigate this effect also for domain walls (or higher dimensional branes). 
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APPENDIX A: DETAILS OF THE CALCULATION OF HIGHER ORDER PERTURBATIONS 

In this Appendix we shall present some of the details of the the 2nd and 3rd order perturbative calculations necessary to find 
the force acting on the kink and its acceleration. 

All equations arising in perturbation theory are 2nd order linear inhomogeneous partial differential equations of the form 

= fW(UW{4> s ),Z®) (Al) 

where the inhomogeneous term must be computed from the solutions obtained in order lower than 0(A n ). To avoid to present 
too complicated general formulae, we start by writing out explicitly the 2nd order equations. They are obtained by simply 
substituting the perturbative expansion for the field cf>(x, t) Eq. il3[ and for the frequency d27l ) into the equation of motion (0: 

P> - 2a/ l V 1 ^ 1 > + = ~U"\4>.) Ue^ ~ 2 W - q vi q e- 2i ^ ■ (A2) 

On the r.h.s. of Eq. (IA2t there are two source terms oscillating with frequency ±2u> and a time independent term. Therefore we 
can seek the solutions of Eq. ( IA2t in the form 

£ (2) (x,t) = $j (x)e 2M + ^1 (x)e~ 2iut + ^ (x, t) . (A3) 

Denoting m-th coefficient in Fourier's expansion of £ of order 0(A n ) by £m\ we obtain the following equations: 

- 2J ( V 1 )£( 1 ) = 0, (A4) 

& 3) +t4 2) = -\u'"mv q ri~ q , (A5) 

(L-4o,(°) 2 )^ = ~U'"{^ ±q . (A6) 

The first equation gives immediately id 1 - 1 ) = 0, i.e. there is no correction to frequency in the first order. Projecting Eq. (IA5t onto 

* "(2) ~ "( < 2) 

the translational mode rjt, and using the identity (L£q \-q t ) = (L??t|Co )* = 0. we obtain 

(^ ) \Vt)=-\(U'"(^)v q V- q h)- (A7) 

As we shall show now the r.h.s. of equation flA7t vanishes precisely for reflectionless potentials. To prove this we shall take the 
derivative of the eigenvalue problem of the linear operator L: 

+ U"(<f> a (x))\ T7,(x) = u 2 r h (x) , (A8) 
then multiply it with rj^ q and integrating leads to 

dx r,_ q (t - lo 2 ) rf q = -(U'"^ s (x))r, q n- q \^ s ) . (A9) 

Integration by parts (over some interval (-L, L)) of the l.h.s. of the above equation and using equation dA8t satisfied by rj_ q 
gives 

(r,' qV '_ q - V-W q )\ L _ L = -(U"'(<l> s )v q V- q \Vt) , (A10) 

where we have also used that the translational mode rj t = <j>' s . Finally it is easy to calculate the boundary values in dAlOb by 
using the asymptotic form for rj q ( fTST l and ( fT~9T > leading to the interesting identity 

4q 2 \n\ 2 = -(U'"(<j )s ) Vq r ] _ q \T lt ). (All) 

This demonstrates that for reflectionless potentials i.e. R = 0, (U'" (4> s )r] q ri- q \r)t) = 0. Therefore it is fully consistent to assume 
^q 2 "* = in Eq. JA51 > when R = 0. It is not difficult to obtain the asymptotic form of £q 2 -* (x) since it has a limit for x — > oo: 

a ^ ±oo) = -4c7^)- (A12) 
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To obtain the asymptotic form of the solution of equation ( IA6t we make use of the Green's function for the inhomogeneous 
problem. Since all inhomogeneous equations we have to solve can be written as 

i^=f^(x), (A13) 
we give their solution satisfying our boundary conditions in a general form: 

= --^^ f dx' Vn{x')f^{x') - ^ Hdx' r,- K (x')f£>&) (A14) 



where r/ K are the solutions of the homogeneous equation (eigenfunctions of L), K = ^/m 2 oj 2 — U"(4> V ac) is the wave number 
corresponding to the frequency muj and W = r} K ri'_ R — if K i]- K = —2in is the Wronskian 2 . To obtain the asymptotic form of the 
solutions we write the integrals as a difference, e.g.: 

dx' ^(aOZ&V) = ( r ~ D dx ' ^0/£V) ■ (M5) 

X. \J — OO J X / 

The integral over the real line can be calculated using the method of residua while the second one can be calculated using the 

in) 

asymptotic form of r/ K and / r „ . This way we obtain the asymptotic form of the solution of Eq. 

A2), _,_ \ U"'((j) vac ) 2 

C + 2{x -> ±00) = nu „^ yi+ q + a 2 2,±k{q)VTk , (A16) 



where k is the wavenumber corresponding to the frequency 2ui, and 

i r°° 

a 2 2.k{q) = - m J dx' mvlV'^s) . (A17) 

((2) (2)' 
£-2 = £+2 

All the coefficients of the homogeneous part will be denoted as 



1 f°° 

W(!)="jy / dx'r, K (x')f^(x'). 



(A18) 



The third order equations have the following forms: 

, L - 9^) = -lu {iv) (4>s)£l 3 ~ U"\^Ml , (M9) 



(A20) 

(n) 



Taking into account the asymptotic forms of £„" equation dA20t can be rewritten as 

L - c(°) 2 ) t[ 3 \x -> ±oo) = ~U'"a 2 2,±k(<l)V- q r)^k- 



5U'" 2 



(A21) 



where derivatives of U has to be taken at <fi = 4> vac - Note that the l.h.s. of the above equation is the same as for harmonic 
oscillator with frequency q. On the r.h.s. there is a source term which oscillates with the resonant frequency. Therefore the 
following condition must be fulfilled to cancel this resonance term 

w(°) \48U" 16 J ' 



2 Note that as there is no first order space derivative in the equations of motion, the Wronskian is independent of x, and thus can be calculated from the 
asymptotic form of r\ q and rj— q . 
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This gives the first correction to the frequency. (Note that the values of the derivatives of the potential must be taken at vacuum.) 
Having this we can write the asymptotic form of the eq. dA20b in much simpler form 

(£ - ^ (0)2 ) £l(x ±00) = -^U"'a 22i±k (q)r)- g ri^ k (A23) 

which leads to the solution 

_i_ ^ ot22,kU"' 

C±'(x -> ±00) = — ^VTkV-q + Ot3-L,±qV^q, (A24) 

p± + U" — U}(°) 

where p± = ^fk — q. A computation similar to the above yields 



d 3) (x ±00) = — ^— ( U"U ( '^ + U'" 2 ) 



1 Cr"'a 22|fc M (A25) 

To — 7772 V^fkVq + a 3 3,± S ??Ts> 

2p|-9w(°) +C/" 

where s is the wavenumber corresponding to the frequency 3a\ As one can see in the third order solution there are only terms 
which oscillate with time, therefore its projection onto the translational mode gives no contribution (after averaging in time) to 
the time independent part of the acceleration. The next time independent term appears at 4th order, 0(A 4 ). Then, the equation 
for ^g 4 "* has the form: 



(2) 2 



(A26) 



Now computing the projection of ^q 4 -* to the translational mode, rj t , exploiting the obvious reflection symmetries some of the 
terms in Eq. jA26i will not contribute, and we are led to the following result: 

(i^\ Vt ) = -Re([/"'(0 s )(^V 9 + i { 2 2) £l) + \u {iv \4> s )£, { 2 2) v 2 - q \m) ■ ( A27 > 



APPENDIX B: NEWTON'S LAW 



In this section we show by explicit computation that Newton's law holds to leading order in PT. We start by calculating the 
projection of the second time derivative of the perturbation £ onto the translational mode of the kink. Computing the second 
order acceleration in PT, by projecting £ on the translational mode, one obtains 



(ma) 2 



L ti(u"\<t> s )d 1 k^)dx=u'\<t> s )d 1) £l L 

L - 1 



u"( (t >M i) 'ei+d i) ei') 



(Bl) 



Next looking at the relevant component of the stress-energy tensor, 



(B2) 



one sees that the first and the last terms give negligible contributions (—Tn)\_ L , if L is sufficiently large. The second term is 
exactly the boundary term in (IB 11 1. To prove the equality, let us put 



-L 



and note that = iuj^ 1 ^. Therefore, the remaining terms in (—Ti{)2\ I L L are equal to 



-L 



^ y u"(<t> s )ei+etu"(^\ 
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where we used the equation of motion L£±i = 0. Comparing this with the last term in (IBU completes the proof. 

If the acceleration is vanishing to second order, the first nontrivial contribution to it can only come form the 4th order. 
This is because neither ma nor Tn has a zero frequency part in the third order. The fourth order acceleration is given by 
{m s a)i = — (Co 4 "*) where £q 4 ' can be replaced by its source term (see equation dA26l )). This should be equal to (— Ti\)i\^_ L , 
where 

(-T U )4 =u"{*.){\g> 2 +g ) e\ +m 3) +^ei} + ^u^pa^ + lei 2 ^ + ^ i)2 a} 

We remark that to remove the resonance terms from the source of £±l, one has to introduce a nonzero . This term induces a 
fourth order correction to (Tn)2. Explicitly it is given by 

To establish Newton's law to this order proceeds essentially the same way as in the second order case. Performing the partial 
integrations in the integral (7710)4 gives all the terms in the energy-momentum tensor containing the potential as boundary terms. 
The integral terms after the partial integration in (7770)4 are 

-/ £ J^(^)(^Wi a )'- 

Note, that the terms in the integrands are the source terms of equations of motion (fT4l . ( IA5b , ( IA6b and (IA20b . Replacing them 
with the right hand side of these equations of motion, and reorganizing the terms, one gets total derivatives, the integrals of 
which are the remaining terms of (— Tn)|£ L . 

Formula ( |B3t can be used to obtain the momentum in the segment [—L,L] to order 0{A 4 ). The rate of change of the 
momentum flow is given by Eq. (fTOb and in fact we need its integrated form Eq. (|28| i. Direct substitution of the series (fT3] l and 
expansion into powers of A shows that for reflectionless potentials the first non-vanishing term is of order 0(A 4 ), therefore 
dtP = A 4 dtP^- We need to calculate the energy-momentum tensor at the boundaries of the segment Tn proportional to A 4 , 
— |£ l . It is sufficient to use the asymptotic form of the solutions ( IA12b . ( IAI6I 1, ( lA24b and JA251 > together with the frequency 
correction (IA22b . The calculation of the boundary term leads to a complicated expression, which after averaging in time gives 
the leading term for the force: 

F = -A 4 [2k 2 (|c4 i+fc | - |al 2> _ fc |) - 4g 2 Rea 3 i,_ (? ] . (B5) 

At the end of this appendix, we will show that to leading order the energy of the accelerating kink is ^mv 2 indeed, as used in 
Section |III] Writing <f> = <f) s + £, and using Eq. (0 we obtain 

e = l/s 2 + urn + + u\<t> e )t + \e + \e + Iv^e + \u"\<t>.)e + ^t/ (iv) (^ 4 ■ w 

The terms linear in £ give no contribution because <j) s solves the static equation of motion. The only time-dependent non periodic 

terms that (up to the fourth order) are those containing £q 4 ' = — \at 2 (j)' s + . . . . In the fourth order, there are no such terms. In 

(2) 

the sixth order, the terms quadratic in £ give a contribution that is proportional to the equation of motion of £q ' multiplied by 
^q 4 "* (in one term one has to perform a partial integration). The term \U'" (</> s )£ 3 gives a term proportional to C^'C-iC^Vsj 
which is antisymmetric, therefore its integral is 0. The last term contains no sixth order contribution proportional to £q 4 ' . Finally 
in the eighth order, the quadratic terms give the contribution \mv 2 . 
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APPENDIX C: SUPERPOSITION OF THE SG KINK WITH A CNOIDAL WAVE 

In this Appendix we review first the remarkable analytic solution corresponding to the nonlinear superposition of a sG kink 
and a traveling wave. This solution has been obtained by H. J. Shin by Darboux transformation methods ifioll . In Ref. ifioll light 
cone coordinates are used, and the sine-Gordon equation is written as: 

d s d z <j> = 2/3 sin 20, (CI) 

where the scale parameter, (3 is kept for bookkeeping purposes as in Ref. ifioll . In these coordinates the static kink takes the form 



X = 2 arctan e - 2 ^( z +*) . (C2) 



The transformation between the two conventions is: 



-x + t _ —x — t 

z 



_ V£ ' (C3) 

X = -2y/0(z + Z) t = 2y/0(z-z). 

and 

<j>(x,t)=24>(z,z). (C4) 
The solution of interest for our purposes is called type 2 in Ref. fiotl . and it can be written as: 

t> - 2k^fl cnfe, *) + l^ ^H ^L. , ,C5) 

and 

sin 20 = £7(1 - fc 2 sn( X2 , k 2 ) 2 ) + Vk sn( X2 , k 2 ) , (C6) 
where T = 1/(5 + 1/S), U = AT {1 - 2T/S) and V = -2(1 - 8T 2 ), and 

X 2 = 2^(z - Vz) (CI) 

and 

^ _ aksn(u,k 2 )cn(x2~u,k 2 )X + bdn(u,k 2 )y 
bk sn(u, k 2 ) cn( X 2 + u, k 2 )y + adn(u, k 2 )X ' 

with 

X = cxp(Mz + fciV X2 )6t(x2 - u), y = exp(-Mz - fciV X2 )e t ( X 2 + u), (C9) 



where 



M 



N 



cn(it, k 2 



ke t (u) ^ 2k dn(w,fc 2 ) tm(u, k 2 ) 



dn(u, k 2 ) sn(tt, k 2 



sn(u, k 2 ) dn(w, k 2 ) 
cn(u, k 2 ) 



2k 



2 sn(u, k 2 ) cn(u, fc 2 )" 1 



cn(u, fc 2 ) 
T Sn ( u > k 2 ) cn(u, fc 2 ) 
dn(w, fc 2 ) 



dn(it, fc 2 ) 



(CIO) 



e t (u) = e 4 



2(K-iK') 



= 1 



2j2(-) n q n ' 



TITTU 
COS( ), 

K K- iK' h 



(Cll) 



with q = cxp[— ttK' / (K — iK')], K = K(k 2 ) and K' = K(l — fc 2 ). Here, sn, cn and dn denotes Jacobi's elliptic functions, 
and K is the complete elliptic integral. We use here the notations and conventions of Abramowitz and Stegun IJJJ]. 

We would like to point out that the parameter k in Eq. ( lC5t should not be confused with the wave number, also denoted by fc 
in the previous sections of this paper. In fact as it turns out k = A/ 2, and therefore it is a suitable expansion parameter. As it 
stands O t (u) is not well suited to be expanded for fc — * 0. It is more convenient to transform it in a form that is easier to expand 



21 



for small k, when q — > for fc — > 0. In fact applying a modular transformation for the function 6*4(2, g), 8t(w) can be brought 
to the form 



t (u) =e l T 



if - iiT' 



K 



■ exp 



4(K-iK')K' 4KK' 



mi 
2K 



iK' 
~~K 



where q = exp(iirr) and it is now straightforward to expand the above function for k — > 0. Let 

iK' 



G(u) 



ITU 

2K 



K 



(C12) 



(C13) 



Since S depends only on X jy, the common singular parts can be dropped. The remaining functions are analytic in k. Thus 



and 



kN' = kN - 



X' = cxp(Mz + kN' X 2)84 
y' = cxp(-Afz - kN' X 2)04 



e'(u) 



tt(x2 - u) iK' 
2K T ~~ ~~K 
tt(X2 + u) _ iK' 
2K T ~~ ~~K 



(C14) 



cn(w, fc 2 ) 



2{K-iK')K' 4KK'J Q(u) 2 dn(u, fc 2 ) sn(u, fc 2 ) 



- fc 



2 sn(u, fc 2 ) cn(u, fc 2 ) 



dn(it, fc 2 ) 



(C15) 



Shin's solution depends on parameters u, fc, V and -. As we shall see in the next section, the perturbative solution depends 
on the amplitude A of the incident wave, the initial position xq = of the kink and its initial velocity Vq = 0. In order to obtain 
the dependence of Shin's parameters on the parameters of the scattering problem, we'll do a series expansion in fc. 



1. Zeroth order 

The zeroth order part of the functions in d(f> is 



M = VW[cotw + tanu] , (kN') = ^ cot u, (C16) 



therefore 



Xq = cxp [ y/Vfitaiiuz + \f cot uz j =: e x 



thus 



X = cxp y-y^Vptunuz - y ^ cotuzj = e~ x , ( C17 ) 
5 = -^ = --exp(-2X) ) 



/3 ^e- 2 * 

4£cotM ° — , (C18) 



V 1 + 4e-« 

a 12 

which can be compared to the derivative of a Lorentz boosted kink at .x 

. X — Xo — vt 

Xv = 4 arctan cxp 

This comparison shows that 



vT^v 2 



l 0g b - = ^=, Wtznu=^±Z=. (C19) 

a VI - v 2 VI - v 2 
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In what follows, we will look at the case in which the kink is at rest in the origin, thus xq = 0, v = and therefore a = b = 1 
and W tan u = 1. Of course, higher order contributions to these formulas are possible. Note, that X = —x/2. 
For the calculation of the zeroth order term in sin 20, one needs 

_ 4S (S 2 - 1) 

therefore 

. . „ ~. TT 2 tanhx 

(sin20) o = [/ o = — , (C21) 

cosnx 

in which we have already plugged in the above parameters. This is in agreement with <p = \ = 2 arctan e x . 

2. First order waveform 

M and kN' do not have first order terms. Therefore, neither do X and y. Thus 

d z (j) = 2k J — - 2 J — i- — - 2fccotM tan ) 1 ^ [e x sin«cos(x 2 - it) - e~ 2X smucos(x 2 +«)] , (C22) 

V v V v cosh 2a cosh 2a l j 



which can be written as 



provided that \i = + i.e. 



Thus 



Here we can see, that the amplitude is 



Similarly, sin 20's first order term is 



V = cj + g -U = u-q. (C24) 



= X + 2fci ( % e^ + jy-ge-^ 4 ) . (C25) 



A = 2k. (C26) 



4(1 - 6S 2 + S£)Si 



(sin20)i = (7i + ^ siiix, C/i = - ^ , 1 7i;" j " 1 , Vb = -2(1 - 8T 2 ), (C27) 

and then, with the right choice of the signs of the square roots, 

- (sin 20) i cosh 2 a; 



■(sin20)i, (C28) 
2y/l- (sin 20) 2 UUS11 " ~ J 

which can be simplified to 



H=^e tot +?7_ g e- iwt . (C29) 



3. Second order waveform 



Here, one has to expand T = Tq + kT\ + k 2 T2- The resulting formulae are 

Tl -(TTW 51 ' 2 ~ OT^ ' (C30) 
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U 2 = — — — °\\ , ,; 2 , 4 — — ^-=-, Vi = 32T Ti. (C33) 



Let Ay/JJVf{u) = 4 v //37V 7 cn(x2, fc 2 )/(sn(w, fc 2 ) dn(u, fc 2 )) denote the coefficient of T in d<j>. Then 

90 2 = 4^(/ 2 T + / T 2 ), (C31) 

where 

/o = cot ti, 

4u — sin Au 

^ 2 = m • 2 • 

16 sin u 

Similarly, 

(sin 20) 2 = U 2 - 2U sin 2 X 2 + Vi sin x, (C32) 

where 

4S* (9 - 145 2 + S$) Sf - 4(1 - 55 2 - 55 4 + 5 6 )5 2 
(l + 5 2 ) 

Again, with the right choice of the signs in the square roots, 

- (sin20) o (sin20) 2 (sin20) 2 

4(1 - (sin20) 2 ) 3 / 2 2(1 - (sin20) 2 ) 1 /2' 

which can be simplified to get 

- Au + sin Au 1 „ / 2 

8 sin 2m cosh x \ cosh x 

(2e x a;coshx + cos(wt — qx) sin(2u) (e 2x sinftut — qx — 2u) + siniujt — qx + 2u))) 

x — k a — ■ 

1 - 6e 2x + e ix 

the first term shows that the kink moves with a constant velocity of 

, 9 Au — sin Au 

v = -k z , 

8 sin 2u 

which can be cancelled by adding a correction to y/V tan u, i.e. 

r— „ , Au — sin4it „„ 

VKtanw=l-fc 2 . (C36) 

8 sin 2u 

The change in the first zeroth order term caused by the fc 2 part of the parameters is of the second order in k. It cancels the above 
constant velocity. 

Because of A = 2k, A 2 = Ak 2 , the second order part of the solution in the usual form is the half of the above 2 , which, by 
introducing q instead of u and u> can be further simplified: 

A2) = >(2) a** , A2)-^t , A2) = e'^ 2 ^ (iq- tank x) e'^+^Hg - tanhaQ 
? ?2 ? " 2 ?0 16coshx(l + (Z 2 ) 16coshx(l + g 2 ) 

1 / tanhx ^ 

+ c r — i o - 2x 

8 coshx \ 1 + q- 2 

which is exactly the same as d48l i, the solution obtained using the perturbative method. 



(C35) 



4. Motion of the kink 



The time-averaged (i.e. zero frequency part) motion of the zero of the kink can be also calculated in an exact form. The 
condition is <fi = 0, thus sin = 0. This can be written as 

Mz + kN' X 2 = 0. (C38) 

This leads to 

v + 1 z M - kN'2y/Vp 



v - 1 ~ z 2<JJkN> 
This means that the velocity of the kink is a constant v, which is a free parameter of Shin's solution. 



(C39) 
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